The Newtonian limit of fourth-order gravity 



Hans - Jiirgen Schmidt 

O 
O 
(N 

. Universitat Potsdam, Institut fiir Mathematik, Am Neuen Palais 10 



X 



D-14469 Potsdam, Germany, E-mail: lijschmi@rz.uni-potsdam.de 

Abstract 



J> ■ The weak-field slow-motion limit of fourth-order gravity will be 

■ discussed. 

o ■ 
\o 
o . 

Let us consider the gravitational theory defined by the Lagrangian 



gn! Lg = {8nG)-^ [R/2 + {aRijR'^ + (3R^)V) . (1) 

^1 G is Newton's constant, / a coupling length and a and jd numerical param- 

eters. Rij and R are the Ricci tensor and its trace. Introducing the matter 



^ \ Lagrangian Lm and varying Lg -|- one obtains the field equation 

Eij + aHij + l3Gij = SnGTij . (2) 

For a = f] = this reduces to General Relativity Theory. The explicit 
expressions Hij and Gij can be found in STELLE (1978). 

In a well-defined sense, the weak-field slow-motion limit of Einstein's 
theory is just Newton's theory, cf. DAUTCOURT (1964). In the following 
we consider the analogous problem for fourth order gravity eqs. (1), (2). 
For the special cases a = (PECHLANER, SEXL(1966), POLIJEVKTOV- 
NIKOLADZE (1967)), a + 2f3 = (HAVAS (1977), JANKIEWICZ (1981)) 
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and q; + 3/3 = (BORZESZKOWSKI, TREDER, YOURGRAU (1978)) this 
has aheady been done in the past. Cf. also ANANDAN (1983), where torsion 
has been taken into account. 

The slow- motion limit can be equivalently described as the limit c — > oo, 
where c is the velocity of light. In this sense we have to take the hmit G — > 
while G ■ c and I remain constants. Then the energy-momentum tensor Tij 
reduces to the rest mass density p: 

Tij^S^S'^p, (3) 

x*^ — t being the time coordinate. The metric can be written as 

ds^^{l- 2(j))de - (1 + 2V^) {dx^ + dy^ + dz^) . (4) 

Now cqs. (3) and (4) will be inserted into cq. (2). In our approach, products 
and time derivatives of and can be neglected, i.e., 

R = A/^iIj - 2A(p , where A/ = f\^^ + f\yy + f\,, . 

Further Roq = -A0, Hqo = -2Ai?oo - Ai? and Goo = -4Ai?, where 1^1. 

Then it holds: The vahdity of the 00-component and of the trace of eq. 

(2), 

Roo - R/2 + aHoo + (3Goo = 87rGp (5) 

and 

- R-A(a + 3(3)AR^S7rGp, (6) 
imply the validity of the full eq. (2). 

Now, let us discuss eqs. (5) and (6) in more details: Eq. (5) reads 

- A0 - R/2 + a(2AA(f) - AR) - A(3AR = 87rGp . (7) 
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Subtracting one half of eq. (6) yields 

- A0 + 2aAA0 +{a + 2(3) AR = AnGp . (8) 

For a + 2/3 — one obtains 

- (1 - 2q;A)A(/) = AnGp (9) 

and then ip = (f) is a solution of cqs. (5), (6). For all other cases the equations 
for and do not decouple immediately, but, to get equations comparable 
with Poisson's equation we apply A to eq. (6) and continue as follows. 

For a + 3/3 = one gets from eq. (8) 

-(1-2q;A)A0 = 477^(1 + 2q;A/3)p. (10) 

The A-operator applied to the source term in eq. (10) is only due to the 
application of A to the trace, the original equations (5), (6) contain only p 
itself. 

For a — one obtains similarily the equation 

- (1 + 12(3 A) A(/) = AnG{l + 16(3 A)p . (11) 

For all other cases - just the cases not yet covered by the literature - the 
elimination of il^ from the system (5), (6) gives rise to a sixth-order equation 

- (1 + 4(q; + 3/5) A) (1 - 2aA)A(j) = AnG (1 + 2(3a + 8(3) A) p . (12) 

Fourth-order gravity is motivated by quantum-gravity considerations and 
therefore, its long-range behaviour should be the same as in Newton's theory. 
Therefore, the signs of the parameters a, (3 should be chosen to guarantee an 
exponentially vanishing and not an oscillating behaviour of the fourth-order 
terms: 

q;>0, q; + 3/3<0. (13) 
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On the other hand, comparing parts of eq. (12) with the Proca equation it 
makes sense to define the masses 

m2 = (2q;)"^/^ and mo = (-4(0; + 3/?))"^/^ . (14) 

Then (13) requires the masses of the spin 2 and spin gravitons to be real. 

Now, inserting a delta source p — mS into eq. (12) one obtains for (f) the 
same result as STELLE (1978), 

= mr'"''' (1 + exp(— mor)/3 — 4exp(— m2r)/3) . (15) 

To obtain the metric completely one has also to calculate ip. It reads 

-0 — rnr~^ (1 — exp(— mor)/3 — 2 exp(— m2r)/3) . (16) 

For finite values mo and m2 these are both bounded functions, also for r ^ 0. 
In the limits a — > (m2 oo) and a + 3/3 — >^ (mo — > oo) the terms with 
mo and m2 in eqs. (15) and (16) simply vanish. For these cases and 
become unbounded as r ^ 0. 

Inserting (15), (16) into the metric (4), the behaviour of the geodesies 
shall be studied. First, for an estimation of the sign of the gravitational force 

we take a test particle at rest and look whether it starts falling towards the 
centre or not. The result is: for mo < 2m2, gravitation is always attractive, 
and for mo > 2m2 it is attractive for large but repelling for small distances. 
The intermediate case mo = 2m2, i.e., 3q; + 8/3 = 0, is already known to be 
a special one from eq. (12). 

Next, let us study the perihelion advance for distorted circle-like orbits. 
Besides the general relativistic perihelion advance (which vanishes in the 
Newtonian limit) we have an additional one of the following behaviour: For 
r — > and r — >^ oo it vanishes and for r ~ 1/mo and r ~ l/m2 it has local 
maxima, i.e., resonances. 
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Finally, it should be mentioned that the gravitational field of an extended 
body can be obtained by integrating eqs. (15), (16). For a spherically sym- 
metric body the far field is also of the type 

mr""*^ (1 + aexp(— mor) + 6exp(— m2r)) , 

and the factors a and b carry information about the mass distribution inside 
the body. 
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